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Effect of fluid-colloid interactions on the mobility of a thermophoretic 
microswimmer in non-ideal fluids 


Dmitry A. Fedosov", Ankush Sengupta 0 , and Gerhard Gompper 0 

Janus colloids propelled by light, e.g., thermophoretic particles, offer promising prospects as artificial microswimmers. However, 
their swimming behavior and its dependence on fluid properties and fluid-colloid interactions remain poorly understood. Here, 
we investigate the behavior of a thermophoretic Janus colloid in its own temperature gradient using numerical simulations. The 
dissipative particle dynamics method with energy conservation is used to investigate the behavior in non-ideal and ideal-gas 
like fluids for different fluid-colloid interactions, boundary conditions, and temperature-controlling strategies. The fluid-colloid 
interactions appear to have a strong effect on the colloid behavior, since they directly affect heat exchange between the colloid 
surface and the fluid. The simulation results show that a reduction of the heat exchange at the fluid-colloid interface leads to an 
enhancement of colloid’s thermophoretic mobility. The colloid behavior is found to be different in non-ideal and ideal fluids, 
suggesting that fluid compressibility plays a significant role. The flow field around the colloid surface is found to be dominated 
by a source-dipole, in agreement with the recent theoretical and simulation predictions. Finally, different temperature-control 
strategies do not appear to have a strong effect on the colloid’s swimming velocity. 


1 Introduction 

The construction of nano- and micro-machines, which can 
move through a fluid environment, is one of the grand chal¬ 
lenges confronting nanoscience todayi"—. Several strategies 
and physical mechanisms have been employed so far to gener¬ 
ate self-propulsion in a fluid—. One approach is biomimetic, 
where the flagellar propulsion of sperm, bacteria, or cilia is 
recreated with synthetic soft materials and actuators. Some 
recent examples include artificial sperm— and artificial cilia—; 
however, in most of these cases, the machines are rather of 
millimeter than of sub-micrometer size. Nano- to micrometer 
length scales are reached by artificial cilia made from micro¬ 
tubules and motorproteins-2, and by magnetic nano- and mi¬ 
croscrews rotated by an external magnetic field—. Another ap¬ 
proach is physico-chemical, where non-equilibrium concen¬ 
tration fields or temperature distributions in the fluid environ¬ 
ment are generated around the swimmer and are employed 
for propulsion, without any movable parts of the swimmer 
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itself. Here, diffusiophoretic— — and thermophoretic 
Janus colloids have been studied most intensively. For dif¬ 
fusiophoretic swimmers, a semi-spherical cap on the colloidal 
Janus particle catalyzes a reaction in the fluid and thereby gen¬ 
erates a spatially inhomogeneous non-equilibrium distribution 
of reaction agents and products. For thermophoresis, a semi- 
spherical cap on a colloidal Janus particle absorbs light from 
an external light source, and is thereby heated, and generates 
a local temperature gradient. An interesting combination of 
the two has also been investigated, where the onset of the 
chemical reaction is triggered by the external light intensity—. 
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Light-controlled microswimmers have the advantage that their 
motion can be controlled easily by a variation of the light in¬ 
tensity. 

A particularly interesting type of thermophoretic mi¬ 
croswimmers has been suggested by Volpe et al.—. This is 
again a colloidal particle with a metallic, light-absorbing cap; 
however, this Janus colloid is immersed in a binary fluid mix¬ 
ture at an ambient temperature just below its lower demixing 
critical point. A slight heating of the cap then leads to a lo¬ 
cal demixing of the fluid mixture, which generates the driving 
force for swimming. An advantage of this mechanism is that it 
works for much smaller power of the light source than that for 
thermophoretic microswimmers in single-component fluids. 

The study of thermophoretic swimming in binary fluids^ 
demonstrates that the interaction of the fluid with the colloid 
surface plays an important role. In particular, it was shown by 
Volpe et al.— that the swimming direction depends on which 
of the two components partially wets the colloid surface. 
However, theoretical studies have only considered either fluids 
on the level of the incompressible Navier-Stokes equation—, 
or fluids with an ideal-gas equation of state 20,22 - 24 . Thus, as 
a first step toward an understanding of thermophoretic self¬ 
propulsion in real fluid mixtures, we investigate a system with 
a single-component fluid with a non-ideal equation of state 
and a variety of boundary conditions (BCs) on the colloid 
surface. In particular, we investigate the influence of no-slip 
and slip BCs and fluid-colloid interactions. In addition, the 
behavior of a thermophoretic swimmer is studied for differ¬ 
ent temperature controls and gradients. We find that local 
fluid-colloid interactions and temperature gradients near the 
colloid’s surface control its swimming velocity. The models 
with non-ideal and ideal fluids lead to qualitatively different 
trends in the colloid mobility. Finally, the flow field around a 






swimming thermophoretic colloid is mainly determined by the 
source-dipole term in agreement with recent theoretical pre¬ 
dictions^-. 

2 Models and methods 

We consider a spherical Janus colloid immersed in a single¬ 
component fluid. In order to correctly describe the hy¬ 
drodynamics of a thermophoretic microswimmer, the fluid 
model has to properly conserve momentum and energy locally. 
We employ here the version of dissipative-particle-dynamics 
(DPD) approach— -— with energy conservation^ 7 ™. 

2.1 Dissipative particle dynamics with energy conserva¬ 
tion (eDPD) 

In the standard (isothermal) DPD approach — the fluid is 
modeled by a set of N particles, each of mass m, interacting 
through a weak conservative (F f), a dissipative (Fy), and ran¬ 
dom (F 77 ) pairwise forces, where the subscripts i. j refer to 
particle indices. The pairwise additive forces are given by 

¥f j = a(0 C {r ij )T i j, 

F/y = - TO) 0 {nj ) (Vy • fy)fy , (1) 

Fy = 

where ry = r,- — r ; , fy = ry/ry is its unit vector, and Vy = 
Vj — V; with r and v being the particle positions and veloci¬ 
ties, respectively. The parameter a is the conservative force 
coefficient, which affects fluid compressibility, y and a are 
the friction and noise amplitudes, which are related through 
the fluctuation-dissipation balance^ as a 2 = 2yksT, where 
ks is the Boltzmann constant and T is the temperature. All 
forces are short-ranged with a cutoff radius r c , and vanish 
for ry > r c . The conservative force profile is defined by 
co c (rjj) = (1 — ry/r c ) for ry < r c , while the spatial depen¬ 
dence of the dissipative (ft) D ) and random (® R ) weight func¬ 
tions is determined by © fi (ry) = (1 — rij/r c ) s for ry < r c 
and the relation co D = (co R ) 2 derived from the fluctuation- 
dissipation theorem—. Here, s is the exponent which affects 
inter-particle friction and fluid viscosity such that s < 1 is ad¬ 
vantageous in order to achieve a sufficiently large fluid vis¬ 
cosity 29,30 . |y is a Gaussian distributed random variable with 
zero mean and unit variance with the requirement <jjy = ty, for 
momentum conservation. Finally, At is the simulation time 
step. 

The particle dynamics is determined by the equations of 
motion as 

f,=vy v^I^FC+F^ + Ff,), (2) 


which are integrated using the velocity-Verlet algorithm—. 

In the energy-conserving DPD (eDPD) metho d - 27 ' —, each 
fluid particle i, in addition to its position and velocity, also 
possesses an internal energy £,. We assume that the internal 
energy of a particle is related to the temperature value 7} as 
£,■ = C v 7], where C v is the specific heat of a fluid—. Then, the 
evolution equation for particle temperature can be written in 
the following form 

C v t i = q i = '£(^ i j + q^), 0 ) 

j 

where q, is the heat flux between particle i and the neighboring 
particles j within the cutoff radius r c . The heat flux q-, is a 
sum of pairwise contributions from the heat conduction (q c -j) 
due to local temperature gradients between particles and the 
work (q^'j) done by the conservative, dissipative, and random 
forces. For instance, the work done by the dissipative force 
corresponds to viscous heating. 

Local heat conduction between particles is defined a s 27 ' 28 

9y = kijCOcirij) ^ ^ + ay®,(ry)CyAr 1/2 , (4) 

where ky is the thermal conductivity coefficient, ay is the 
noise amplitude, and £y is the associated noise modeled from 
the Gaussian probability distribution with zero mean and unit 
variance and with the condition = — Qj. Note that q l \: in 
Eq. 0 consists of deterministic and random heat-conduction 
terms and provides local energy conservation. The conductiv¬ 
ity coefficients are defined as2742 
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k U = (Ti + Tj ) 2 , « 2 - = Zkskij, ( 5 ) 

where ko is a nominal strength of interparticle heat conductiv¬ 
ity. For simplicity, we also select ft) c (ry) = o)^(ry). 

To connect particle dynamics to its internal temperature 
(or energy) 27 ' 28 , the force coefficients in Eq. 0 become 
functions of temperature and should be replaced by the cor¬ 
responding dij(Ti,Tj), jij{Ti,Tj), and Oij(Tj,Tj) coefficients. 
Recently, it has been suggested that the conservative force 
coefficient ay should depend linearly on temperature in or¬ 
der to properly reproduce fluid compressibility—. However, 
in the current work we employ a constant conservative force 
coefficient such that ay(7/,r ; ) = a /•-. For simplicity, we also 
assume no temperature dependence of the random force co¬ 
efficient such that Gij(Tj,Tj) = a. Then, the corresponding 
friction coefficient is given by 

(6) 

The above expression for the inter-particle friction also deter¬ 
mines the fluid’s viscosity, which is controlled by setting the 
parameter o in simulations. 





The work fluxes q'^ can be derived from the total energy 
£ of a simulated system. The total energy of the system, 
E = £,£,■ +E mec h, should be constant and thus, a change in 
mechanical energy E mec h should correspond to the change in 
the internal energy such that dE, nec h = -d (J_ i £,). The me¬ 
chanical energy consists of kinetic and potential contributions 
given by 



where 0(r, ; -) is the interaction potential giving rise to the 
conservative force F fAnj) = ^ 7< l ) (, r ij )■ Therefore, we ob¬ 
tain 0(r) = 0.5ar c (l — r/r c ) 2 . One way 27,28 to derive the 
work fluxes qfj is to take the differential dE mec f, from the me¬ 
chanical energy in Eq. 0 and substitute the corresponding 
terms with the force expressions from Eq. 0 using the equa¬ 
tions of motion. Such a method does not strictly provide en¬ 
ergy conservation, and therefore, relatively small time steps 
might be required to have acceptably small energy variations. 
Another method is to implement energy conservation locally 
without explicit calculation of the q"j fluxes. Here, we as¬ 
sume that a change in the internal energy locally is equal to 
the change in both kinetic and potential energies such that 
dEi = —d(mvj)/2 — d (L/0( r iy)/2) . Thus, we calculate the 
changes in kinetic and potential energies locally after each in¬ 
tegration step and adjust the internal energy of each particle 
accordingly. This method leads to exact conservation of the 
total energy of a system, and we have verified that it properly 
represents temperature gradients in a fluid. A similar idea has 
been used in Ref. 3 -, where energy conservation has been im¬ 
plemented locally for each pair of interacting fluid particles. 

2.2 Janus colloid and boundary conditions 

The Janus colloid is modeled with N s = 4000 DPD particles 
placed on the surface of a sphere with radius R s = 4r c corre¬ 
sponding to the surface density of p s « 20/r 2 . The colloid is 
centered at the origin, with the hot (cold) side positioned in 
the half-space x < 0 (x > 0). The colloid particles are frozen 
at their positions on the spherical surface, see Fig.|T| It is im¬ 
portant to note that a colloid moving in a resting fluid and a 
colloid fixed at a position with certain orientation are essen¬ 
tially identical; the main difference is just a transformation of 
the reference frame. In the former case, the colloid is mov¬ 
ing in the fluid, while in the latter case the fluid is moving 
around the colloid. The two cases are equivalent for small 
rotational and translational diffusion coefficients, i.e. for suf¬ 
ficiently large colloids. There might be deviations for small 
colloids, when the relaxation time of the temperature profile 
becomes comparable to the rotational diffusion time. 

In the simulations, one or both sides of the colloid surface 
are maintained at constant temperature. The steady-state be- 



Fig. 1: Schematic diagram of the Janus colloid with the A and B 
caps. The inset is a zoomed part of the cap along an arc showing 
two types of immobile particles which constitute the colloid surface. 
The cap particles (type A and type B) shown in different colors inter¬ 
act with nearby fluid particles (e.g. i) through the soft short-ranged 
DPD interactions, which are similar to those between fluid particles 
(shown in blue). In addition to the DPD interactions, bounce-back 
and specular reflection boundary conditions are used (here shown in 
the same picture, for convenience) at the colloid-fluid boundary rep¬ 
resented by the arc. For specular reflection, a fluid particle, e.g. j, 
will be reflected (shown by the arrows) at the spherical surface such 
that the tangential velocity component is unaffected and the normal 
component is reversed. For the case of bounce-back condition, a fluid 
particle, e.g. k , will be reflected back along its incident direction as 
shown by the arrows. 

havior of the system is ensured by the constancy of the tem¬ 
perature gradient in the fluid and of the flow field in the co¬ 
moving frame of the colloid. The fluid particles are present 
both outside and inside the colloidal shell, in order to ensure 
a proper pressure balance on both sides of the colloid surface. 
Thus, the fluid particles inside and outside the colloid interact 
through the conservative force in Eq. 0. However, the fric¬ 
tion coefficient between the inside and outside fluids is set to 
zero, since viscous interactions between them are shielded by 
the colloid wall. This means that jij = Gij = 0 whenever (i.j) 
represents a pair of fluid particles intercepted by the colloid 
surface. 

Two separate cases of boundary conditions (BCs) are con¬ 
sidered including either bounce-back collisions or specular re¬ 
flection of the fluid particles at the fluid-colloid interface. This 
applies to collisions both at the inside and the outside of the 
colloidal shell. The reflections of fluid particles are necessary 
to prevent the entry of a fluid particle from the exterior to the 
interior of the colloidal shell, or vice versa, since the conser- 






vative interactions between fluid and colloid particles are too 
soft to guarantee no inter-fluid mixing. On the other hand, 
the two different collision rules realize different BCs at the 
colloid surface. The bounce-back collisions implement a no¬ 
slip (or stick) BCs, since the particle velocity is inverted at the 
colloid surface (i.e. v — > — v) resulting on average in a van¬ 
ishing tangential component of the fluid velocity. In the case 
of specular reflections, the velocity component parallel to the 
local tangent plane of the colloid surface remains unchanged, 
while the perpendicular component is inverted. This realizes 
slip BCs. 

In addition to the above mentioned BCs, the fluid parti¬ 
cles also interact pairwise with the immobile DPD particles 
comprising the colloid surface. Therefore, different halves of 
the Janus particle possess not only dissimilar thermal proper¬ 
ties, but also may have different fluid-colloid interactions (see 
Fig.ffll. This is done by assigning different pair-interaction 
coefficients for the particles at the two halves of the colloid 
with the surrounding fluid particles. The particles of the two 
distinct hemispheres are labeled as A and B. Given that the 
fluid particles i and j interact with a coefficient a, ; - = a/,-, the 
coefficients for fluid-colloid interactions will be referred to as 
ciiA and cijB for the two halves, respectively. Then, a u \ < ay 
mimics an effectively “solvophilic” surface, while a u \ > ay 
mimics an effectively “solvophobic” surface. In experiments, 
Janus particles can be chemically functionalized to generate 
hydrophilic or hydrophobic interaction with the host fluid—. 
Such chemical functionalization can affect the particle behav¬ 
ior in addition to the temperature gradient. Thus, changing 
the fluid-cap interaction strength for both caps A and B al¬ 
lows us to access surface-tuning capabilities together with the 
thermophoretic control in our model. The dissipative and ran¬ 
dom force coefficients for fluid-colloid interactions are set to 
zero (jiA = <7,..I = YiB = <7,g = 0), since the bounce-back and 
specular reflections of particles at the colloid interface already 
define the type of BCs employed. Finally, the fluid and the 
colloid exchange heat locally within r c following Eq. (0). In 
order that the fluid particles can approach the colloid particles 
close enough and exchange heat, a shorter cut-off radius for 
the fluid-colloid conservative interactions, r' c = 0.25 r c , has 
been used. 

The temperature gradient across the colloid is maintained 
by setting an elevated temperature Tj wt for the immobile par¬ 
ticles comprising cap A, while the particles comprising cap B 
are kept at a lower temperature T co m throughout the course 
of the simulation, see Fig. Eta). In experiments, however, 
one-half of the Janus colloid is heated by allowing the metal- 
capped half to absorb heat from an incident laser fiel d 17 ' 35 . 
The bulk fluid in this case remains at a colder temperature, and 
thus maintains a temperature gradient across the colloid diam¬ 
eter without the need to additionally cool the other half of the 
colloid. Therefore, we have also examined the case where we 


keep the particles of cap A at an elevated temperature 7} )of and 
impose the temperature T cn i r j for the fluid particles appearing 
within a narrow slice of 1.5r c at the two periodic boundaries 
parallel to the asymmetry axis of the Janus colloid (Fig. [2b)). 
In this case, we let the temperature of all other particles, in¬ 
cluding those of cap B, relax to an intermediate steady-state 
value. Similarly, we have also compared these cases with the 
situation when all the periodic boundaries are kept at the lower 
temperature T co id within a slice of 1.5 r c , while the particles of 
cap A are maintained at the higher temperature Ti wt (Fig.EJc)). 
We denote a maximal temperature difference across the col¬ 
loid as AT = T,\ 7/;, which is equal to 7) 10/ — T co i ( j for the 

case where the temperature is controlled on both sides of the 
Janus colloid. For the other two cases, where the temperature 
is maintained only at one side of the colloid, AT = T] lot — 7g 
with 7g measured directly from the simulation data. The dif¬ 
ferences in Janus-colloid behavior with respect to the different 
temperature-control strategies will be discussed. 

2.3 Simulation setup and parameters 

In the simulations, a cubic box of size L x = L y = L z = 25r c 
with periodic BCs is used. We employ m = 1 to define 
the unit of mass, r c = 1 to define the unit of length, and 
T = r c yjm/{ksT) = 1.25 is the unit of time, where kg = 1 
and T = (Tf, ot + T co id)/2 is the average temperature, which 
has been set to ? = 0.64 in simulations. For the normal¬ 
ization of different colloid properties, we will also use the 
rotational diffusion coefficient D, = kgT /(8tn] R ]), where 1 ] 
is the fluid’s dynamic viscosity at T. Other parameters in¬ 
clude the density of fluid particles Pf = 3/r 3 , the conserva¬ 
tive force coefficient a/.- = 39kg? /r c , the random force coef¬ 
ficient Cjj = 4.2kg? \/T/r c , the exponent s = 0.25, the spe¬ 
cific heat C v = 200kg, and the nominal strength of interparti¬ 
cle heat conductivity ko = 0.00125/t. The fluid viscosity for 
dF = 39kg? /r c at ? is equal to rj = 2.74i/mkg? /r 7 , and thus 
D r = 2.27 x 10 -4 i/kg ?/(mr 7 ). We will also employ a fluid 
with dp = 0 whose viscosity is equal to rj = 1.78 y/mkg? /r 7 
resulting in D, = 3.5 x 10~ 4 y/kg? /(mr~). The simulations 
were performed with discrete timesteps Af = 0.004 t and data 
were collected after an initial relaxation time of 2 x 10 s 
timesteps, in order to ensure a proper steady state. The re¬ 
laxation to the steady state has been monitored by observing 
a stable time-independent temperature gradient and flow field. 
The steady-state averages were calculated by collecting data 
over another 8 x 10 6 timesteps. 

3 Results 

We study the dynamics of the model thermophoretic Janus col¬ 
loid in the host fluid medium under various temperature condi- 
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Fig. 2: Steady-state temperature profiles within the two-dimensional 
cross-section of a spherical Janus colloid (black semi-circle) for 
(a) a differentially heated surface of the colloid with temperatures 
Thot/T = 1.3 and T co id/T = 0.7 maintained on the hot and cold caps 
of the colloid, respectively, (b) the colloid heated only at one hemi¬ 
spherical cap maintained at 7}, of /? =1.3 and with a temperature 
control at the periodic boundaries in the y—direction maintained at 
Tcoia/T = 0.7, and (c) the colloid heated only at one hemispheri¬ 
cal cap maintained at Ti w ,/T = 1.3 and with the temperature con¬ 
trol at all periodic boundaries maintained at T co u/f = 0.7. Here, 
T = (Ti lot 4- T co id)/2 is the average temperature, which has been set 
to T = 0.64 in simulations. These simulations correspond to the case 
of af = 39 kgf /r c and cijA = UiB = 7.8 hgt/r c . 


tions and surface interactions. The center of the Janus colloid 
is fixed at the origin of the reference frame as discussed in 
Sec. l2.2l such that we can conveniently study the translational 
motion of the colloid by measuring the fluid velocity far from 
the colloid surface. All our results are shown for steady-state 
conditions. 


3.1 Temperature dependence of self-propulsion 

Figure [2] shows the temperature profiles around the colloid in 
the steady-state for the different temperature-control strategies 
described above. The temperature profiles are axisymmetric 
allowing us to average simulation data over the full azimuthal 
angle. Clearly, the temperature profiles are very different in 
the three cases. In Fig. 0 a), a very strong temperature gradi¬ 
ent develops at the interface between the two caps; this inter¬ 
facial gradient is much smaller in Figs. [2b) and [2c), which 
is both due to the smaller temperature difference between the 
two caps and the temperature variation on the non-heated cap 
itself. The comparison of the profiles of Figs. [2b) and [2c) 
shows minor differences in the close vicinity of the colloid 
surface, but a much slower decay along the symmetry axis in 

Fig. 0b). 

In all three cases, we also measure the density profiles of 
the fluid near the colloid surface. Figure [2 shows the fluid 
density cuts along the x— axis and the corresponding density 
profiles in two dimensions. We find a layering of the interact¬ 
ing fluid particles at the fluid-colloid surface. Such a layer¬ 
ing is well known for hard-core particles near a hard wall—, 
and is thus related to the repulsive interaction of fluid parti¬ 
cles among themselves and with the colloid surface. These 
density modulations at the colloid wall decay to the mean 
bulk density (py = 3 /r%) rapidly over a length scale of two 
to three times r c . However, the density patterns are similar 
for different temperature-control strategies. Also, density pro¬ 
files do not show any appreciable differences for the various 
BCs, in particular the bounce-back and the specular reflection 
BCs at the fluid-colloid interface. The fluid density profiles in 
case of a hot cap (7} )of /f = 1.3) and a temperature fixed to be 
Tcoid/T = 0.7 at the periodic boundaries in all directions are 
nearly identical to those in Figs. 0b) and0d) for the temper¬ 
ature control only at the periodic boundaries in the y-direction. 

The non-zero temperature difference across the colloid di¬ 
ameter results in the spontaneous generation of a far-field flow 
velocity in the fluid in the co-moving frame of the colloid 
(see Fig. 0). We measure this flow velocity by averaging 
over the velocities of all the DPD fluid particles far away 
from the colloid with coordinates |r| > 2 R s . In the labora¬ 
tory frame, this flow velocity is the same in magnitude and 
opposite in direction to the self-propulsion velocity v p of the 
colloid. We find that v p increases linearly with increase in 
the maximal temperature difference AT across the colloid at 
fixed average temperature T, as shown in Fig. 0 in terms of a 
particle Peclet number Pe p = v p /(2R s D r ). Pe p can be also 
interpreted as a non-dimensional swimming velocity of the 
colloid. Defining the corresponding proportionality constant 
as the thermophoretic mobility, jx = v p /AT, we can readily 
obtain /I from the slope. The thermophoretic mobility is de¬ 
termined by the surface properties of the Janus particle, the 
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Fig. 3: Density (p) profile of the fluid around the Janus colloid in the steady-state. The variation of the fluid density extracted along the x —axis 
at r = 0 is shown for the Janus colloid with (a) a hot = 1.3) and a cold (J co id/T = 0.7) cap and (b) a hot cap (Tp ot /T = 1.3) and 

a temperature fixed to be T co u/T = 0.7 at the periodic boundaries in the y— direction, (c) and (d) show the corresponding density profiles 
in two dimensions averaged axisymmetrically over the full azimuthal angle. These simulations correspond to the case of ap = 39 kpT/r c , 
cijA = cijB = 7.8 kgf /r c , and bounce-back reflection BCs. 


interactions within the fluid, and the average temperature T, 
but obviously independent of AT, and thus, it is a convenient 
quantity to characterize a thermophoretic microswimmer. We 
can then use the thermophoretic mobility fi to study the depen¬ 
dence of self-propulsion on the surface properties of the col¬ 
loid in relation to the host fluid. Subsequently, our simulation 
results will be presented mostly in terms of a non-dimensional 
mobility /J* = Pe p T/AT = /if/(2 R s D r ). 

3.2 Dependence of thermophoretic mobility on fluid- 
colloid interactions 

3.2.1 Effect of fluid-colloid repulsion. The surface 
properties of the Janus colloid can be tuned to manipulate and 
control its thermophoretic mobility. We first consider the case 
of “symmetric caps”, when — apart from the different tem¬ 
perature conditions at the two caps (T& > Tb) — the surface 
interaction of both caps A and B with the fluid particles i are 
identical, i.e., a, a = a iH = Ucap. The dependence of the self¬ 
propulsion velocity on AT is shown in Fig. [fija) for a wide 
range of cap interactions, with 0 < a ca p/aF < 8. These results 
indicate that the thermophoretic mobility depends highly non- 
linearly on the fluid-cap repulsion strength, a cap , as shown in 
Fig. [6] for the case of AT/T = 0.6. The mobility increases 
with increasing cap repulsion strength, and levels off to some 



Fig. 4: Streamlines for fluid flow around the thermophoretic 
microswimmer with symmetric cap interaction (a ,a = a iB = 
312 kpT/r c ) in the host fluid (ap = 39 kgf /r c ), and bounce-back BCs 
at the fluid-colloid interface (black semi-circle). The color-code cor¬ 
responds to the v-component of the fluid velocity (v x ) shown in terms 
of the local Peclet number Pe x = v x / (2 R s D r ). This simulation corre¬ 
sponds to the case when the temperature is controlled at the colloid 
surface with a hot ( Tp ot /T = 1.3) and a cold (J co id/T = 0.7) cap. 


saturation value at very high repulsion strength. 

Figure |6]reveals an interesting dependence of jx* on the slip 
or stick BCs. Within the numerical accuracy, we find that the 
saturation value at high repulsion strength does not depend on 


















































(b) a ii =7.8k |J T/r (; in all cases 



a iB =7.8k B T/r c in all cases 



Fig. 5: Self-propulsion velocity v p presented in terms of the particle 
Peclet number Pe p = v p /(2R s D r ) as a function of the temperature 
difference, AT = — Tb, for a fluid with interaction strength cif = 

39 ksf /r c , when (a) a = am (symmetric cap repulsion strengths), 
(b) a,A < ap (an effectively solvophilic hot cap), and (c) am < ap 
(an effectively solvophilic cold cap). All data points are obtained for 
the bounce-back BCs at the fluid-colloid interface. All simulations 
correspond to the case when the temperature is controlled directly at 
both caps of the colloid. 


whether bounce-back or specular-reflection collisions are em¬ 
ployed at the colloid surface. This can easily be understood by 
considering the fact that at high repulsion strengths hardly any 
fluid particles can reach the colloid surface anymore, so that 
the type of surface reflection becomes irrelevant. However, 
for low cap repulsion strengths, the stick BCs generate larger 
fluid flow (and thus, propulsion velocity) than the slip BCs for 
a given temperature difference and cap repulsion strength. In 
fact, in the absence of cap repulsion, the thermophoretic mo¬ 
bility disappears in the case of slip BCs. 

Figure 0 also shows the velocity-temperature graphs for 
asymmetric cap repulsion strengths. Figure0b) concerns the 
case when the hot cap A interacts with a lower interaction co¬ 
efficient with the fluid particles compared to the fluid-fluid in¬ 
teraction (aiA < ap), to mimic a solvophilic cap interaction. 
This interaction is kept fixed and the interaction of the cold 
cap B with the fluid is increased from solvophilic (a,/; < ap) 
to solvophobic (am > «/■) strengths. In Fig. 0c), the cold cap 
B is instead maintained at a solvophilic interaction strength 
(aiB < up), and the interaction a L \ on the hot cap A is var¬ 
ied. These two cases are not identical, because the A-cap is 
always the hot and the B-cap is always the cold side. The ve¬ 
locity response to the temperature difference remains linear 
for all cases. The non-dimensional thermophoretic mobility 
jJL * extracted from these data for the case of AT/T = 0.6 is 
shown in Fig. 0 As for the symmetric case, j.1* is found to in¬ 
crease with increasing repulsion strength a ^ or am. Figs. 0a) 
and 0b) also provide the comparison of mobilities for stick 
and slip BCs, which show similar trends as for the symmetric 
case, with an increase in ji* and its final saturation as the cap 
repulsion strength is elevated. 

The thermophoretic mobility obtained from these two very 
different asymmetric fluid-cap interaction cases remains es¬ 
sentially identical (within our numerical accuracy), although 
it is much lower than the corresponding mobility for symmet¬ 
ric fluid-cap interactions (where the larger value of the two 
repulsion strengths in the asymmetric case is the same as a cap 
in the symmetric case). In particular, the saturation value of 
jit* in the asymmetric case is considerably smaller than that 
for the symmetric case. Furthermore, Fig. 0 shows that the 
increase of the mobility due to cap repulsion is about twice 
as large for the symmetric than for the (highly) asymmetric 
cases. Together, this indicates that the thermophoretic mo¬ 
bility can be understood as a superposition of the mobilities 
generated by the two caps independently. This indicates that 
the temperature gradient between the cap and the fluid, rather 
than that at the interface between the two caps, is responsible 
for the propulsion. Figure 0 also shows that the insensitivity 
of the thermophoretic mobility to the exchange of repulsion 
strengths is also not affected by the type of BCs at the colloid- 
fluid interface. However, the slip BCs always generate lower 
mobilities compared to stick boundary BCs, in particular for 
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Fig. 6: Non-dimensional thermophoretic mobility fi* as a function of fluid-cap repulsion strength with (a) bounce-back and (b) specular 
reflection BCs and for interacting fluid particles with of = 39kg? /r c . The hot and the cold caps (A and B) are taken to be interacting with the 
fluid particles symmetrically with a, A = ci,b (red triangles), and asymmetrically with a, A < a/g f° r aiA = ag/5 and varying cz,g (solvophilic hot 
cap, green squares) and with a > «,g for a,g = ap/5 and varying cij A (solvophilic cold cap, blue circles). All simulations correspond to the 
case when the temperature is controlled at the colloid surface with a hot ( T\ lat jT = 1.3) and a cold (T co id/T = 0.7) cap. 


lower values of fluid-cap interaction strengths. 

3.2.2 Effect of temperature control. Figure [6] reflects 
the dependence of the thermophoretic mobility when both 
sides of the Janus colloid are maintained at fixed tempera¬ 
tures T a = T] wt = 1.3? and 7/ ; = T co m = 0.7? (illustrated in 
Fig. |2j a)). We now consider the case where the A-cap is held 
at a fixed temperature T /,,,,, and the fluid far away from the col¬ 
loid (and near the periodic boundaries of the simulation box) 
is held at a lower value T ro hh while the temperature of the 
other cap ( B) is free to adjust, as illustrated in Figs. [2}b) and 
[2c). In this case, 7g is measured directly from the simulation 
data. The temperature distribution is nearly independent of the 
fluid-colloid interactions and the type of BCs, and the mea¬ 
sured values of temperature at the cold side are 7g = 1.11? 
when the two periodic boundaries parallel to the symmetry 
axis of the Janus colloid are kept at the temperature T l:o i,i and 
Tb = 1.09? when all periodic boundaries are maintained at 
T co id- The results for the non-dimensional thermophoretic mo¬ 
bility displayed in Fig. [7] show that the qualitative trends re¬ 
main unaltered. However, fi* is slightly lower than for the 
case of fixed temperature Tb = T co u. This is due to the fact 
that temperature gradients at the interface between the cold 
and hot sides of the colloid are stronger for the case of the 
temperature control Tb = T co id directly at the colloid. 

3.2.3 Comparison with previous simulations. Over 
the entire range of temperature and interaction potentials con¬ 
sidered, we do not observe any reversal of the propulsion di¬ 
rection. The colloid is always propelled in the direction of the 


non-ideal fluid 



Fig. 7: Non-dimensional thermophoretic mobility p* as a function 
of fluid-cap repulsion strength with bounce-back BCs, for interacting 
fluid particles with ap = 39fcg? /r c and with a temperature gradi¬ 
ent maintained between one of the caps (A) kept at an elevated tem¬ 
perature (Ti wl = 1.3?) and the periodic boundaries maintained at a 
fixed lower temperature (7j, 0 /^ = 0.7?). The green squares represent 
the thermophoretic mobilities when two opposite periodic boundaries 
parallel to the symmetry axis of the Janus colloid are kept at the tem¬ 
perature T co ij, and the blue circles correspond to the case when the 
temperature at all six periodic boundaries is maintained at T co ij. For 
comparison, red triangles represent the case where the temperature 
T co id is maintained directly at cap B. The cap repulsions are taken to 
be symmetric {a, A = a,-#) for all cases. 






















cold side. Recent simulations 20 - 22 ' 24 using another mesoscale 
hydrodynamics simulation technique, the multiparticle colli¬ 
sion dynamics (MPC) metho d 37 ' 38 , have shown that the swim¬ 
ming direction of a thermophoretic colloid can be affected by 
the solvent-colloid interaction potential. In these studie s 20 ’ 24 
only a symmetric case was considered, and an attractive poten¬ 
tial between solvent particles and the thermophoretic colloid 
has led to swimming in the direction of the cold side, while 
a repulsive potential has triggered the colloid to move in the 
direction of the hot side. In our simulations, the conserva¬ 
tive potential is purely repulsive; however, for the conditions 
UiA < ai- and a,/j < a/.-, fluid particles should be effectively at¬ 
tracted to the colloid surface. The current simulation results 
do not lead to a similar behavior of the Janus colloid where 
the swimming direction can be interchanged. It is important 
to note that the two fluid models are not equivalent. The main 
difference between eDPD and MPC fluids is that the MPC 
fluid has the equation of state of an ideal gas, while the eDPD 
fluid is much less compressible due to the presence of conser¬ 
vative interactions if ap > 0. Furthermore, the heat capacity 
and thermal conductivity in these models are not equivalent, 
which will be discussed further in text. 

3.3 Thermophoretic propulsion in “ideal” fluids 

The coefficient ap of the conservative interaction between 
fluid particles plays an important role in determining the com¬ 
pressibility of the fluid, but also affects the diffusion coeffi¬ 
cient and the fluid viscosity. Here, decreasing ap increases 
the fluid compressibility. Further, we investigate the extreme 
case when ap = 0, i.e., the case obtained when the fluid is 
most compressible. In this case, the equation of state is that 
of an ideal gas, similar to the MPC fluid model in Refs. 20 ' 24 . 
The thermophoretic mobility of the Janus colloid is then again 
measured in this host fluid medium for both bounce-back 
and specular reflection BCs, and its dependence on the fluid- 
colloid interaction strength for both symmetric and asymmet¬ 
ric cap cases is determined. 

In the case of bounce-back BCs at the colloid-fluid inter¬ 
face, the thermophoretic mobilities as a function of the fluid- 
cap interaction strength are shown in Fig. [Ha). We find an 
initial decrease of the mobility with increasing cap-fluid in¬ 
teraction strength, followed by a subsequent weak recovery, 
and finally saturation with further increase of the interaction 
strength, for both symmetric and asymmetric fluid-cap inter¬ 
actions. This dependence of the thermophoretic mobility is 
qualitatively different from the behavior of a Janus colloid in 
a non-ideal fluid (with ap > 0). Also, the propulsion veloc¬ 
ity in the ideal fluid for the symmetric cap case remains lower 
than that for asymmetric cap, again in contrast to the non-ideal 
fluid case. Our results show that in addition to an enhancement 
of mobility with cap asymmetry, the magnitude can be further 


controlled by interchanging the temperature of the caps, un¬ 
like the non-ideal fluid scenario. 

For specular-reflection BCs in Fig. [8} hj, we obtain a re¬ 
versal of propulsion direction of the colloid in the ideal fluid 
by tuning the fluid-cap interaction. In the low fluid-cap in¬ 
teraction regime, the computed mobility is negative, which is 
equivalent to swimming in the direction towards the colloid’s 
hot side. With increasing fluid-cap interaction strength, the 
corresponding mobility changes sign and becomes positive, 
i.e. the propulsion direction is inverted and the colloid swims 
toward its cold side. Over the entire range of the fluid-cap 
interactions studied in Fig. [Hb), the mobility increased grad¬ 
ually with the interaction strength before leveling off at very 
large fluid-cap interactions. 

A comparison of Figs. [8ja) and[8]b) reveals that the trends 
of the mobilities over the entire range of cap repulsions for 
specular reflections and bounce-back collisions are quite dif¬ 
ferent. In particular, the saturation value of the mobility for 
the colloid with symmetric caps remains higher than that with 
asymmetric caps for specular reflections unlike the bounce- 
back situation. Furthermore, interchanging the temperatures 
of the asymmetric caps leads to a change in the propulsion 
velocity, quite unlike the situation of non-ideal fluids. 

In comparison to the MPC simulations in Refs. 20 ' 24 , we also 
observe a reversal in the swimming direction by tuning the 
fluid-cap interaction strength. However, the direction reversal 
occurs only for the case of specular reflections (i.e., slip BCs) 
and relatively small fluid-cap interaction strength. In these 
simulations, the fluid-cap interaction can be considered purely 
repulsive, since ap = 0 of the suspending fluid. Thus, in case 
of bounce-back reflections the colloid always swims toward 
the cold side similar to the results for non-ideal fluid case, 
while in case of slip BCs, the colloid moves toward the cold 
side for strong fluid-cap repulsion and toward the hot side for 
weak fluid-cap repulsive interactions. This trend seems to be 
rather opposite to that in Refs. 20 ' 24 , where fluid-colloid repul¬ 
sive interactions lead to a swimming direction toward the hot 
side, while fluid-colloid attractive interactions result in swim¬ 
ming toward the cold side. In addition, the swimming direc¬ 
tion of a Janus thermophoretic colloid has not been affected 
by the type of BCs (i.e., slip or no-slip) in the MPC simula- 


3.4 Flow field around thermophoretic Janus colloids 

As the colloid swims, the flow field generated in the fluid 
can be found from our simulations. Figure [9] shows the flow 
field in the lab frame, obtained by adding the velocity of self¬ 
propulsion v p to the fluid velocity v, since its magnitude has 
an opposite sign to the far-field fluid velocity (see Fig. [4] for 
example). The flow field indicates that the source-dipole con¬ 
tribution dominates, which is consistent with the conclusions 
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(a lA . a iB )k B T/r c (a iA , a iB )k B T/r c 

Fig. 8: Non-dimensional thermophoretic mobility (jl*) as a function of the fluid-cap repulsion strength with (a) bounce-back and (b) specular 
reflection BCs for ideal-gas fluid with ap = 0. The hot and the cold caps (A and B) assume symmetric interactions a ,a = am with the fluid (red 
triangles) as well as asymmetric interactions with a fixed = l.Skgf /r c and am varied (green squares) and with a fixed am = 7.8 kgf/r c 
and aiA varied (blue circles). All simulations correspond to the case when the temperature is controlled at the colloid surface with a hot 
(Jhot/f = 1-3) and a cold (T coId /T = 0.7) cap. 



Fig. 9: Flow field v + v p in the lab frame of a thermophoretic 
microswimmer with symmetric cap interaction (a^ = am = 
312 kgt/r c ) in the host fluid with ap = 39 kgT/r c , and bounce-back 
BCs at the fluid-colloid interface (black semi-circle). The color-code 
corresponds to the x-component of the velocity flow-field shown in 
terms of the local Peclet number Pei = (v x + v p )/(2R s D r ). Flere, the 
temperature T co i d is controlled directly at the colloid surface. 


from theory-- and simulation s;'—— - for a thermophoretic col¬ 
loidal swimmer. The flow field in Fig. [9] is very similar to 
the theoretical predictions for a thin-cap limit— (i.e., the ther¬ 
mal conductivity of the cap does not play a role), supporting 
the validity of the simulation results. The corresponding flow- 
field in MPC simulation s ' ?-— is very close to a source-dipole 
approximation in the theory? ■. 

Finally, we determine the radial component v r h of the ve¬ 
locity field v + v p in spherical coordinates as a function of the 
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Fig. 10: Radial velocity v Ks h of the flow field v + Vp in spherical 
coordinates presented in terms of the local Peclet number Pe ri h = 
v rsph /(2 R s D r ) as a function of distance r sp i, from the center of the col¬ 
loid, obtained along fixed polar angles, 0, with respect to the symme¬ 
try axis of the Janus colloid. The filled and the open symbols are for 
bounce-back and specular reflection BCs, respectively. Circles and 
triangles correspond to 6 = n /4 and 6 = 3tt/ 4 cases, respectively, 
for a DPD fluid with ap = 39kg f / r c , and the diamonds correspond 
to an ideal gas fluid with ap = 0 and 6 = k/ 4. The thick solid line 
shows the power-law behavior ~ r~ ph . The simulations correspond 
to the case of a ^ = am = 312 kgT/r c with the temperature T co i d con¬ 
trolled directly at the colloid surface. 














































Fig. 11: Radial velocity v r h of the flow field v + v p in spherical 
coordinates presented in terms of the local Peclet number Pe r h = 
v rs h /(2R s D r ) as a function of the angle (0) with respect to the sym¬ 
metry axis of the Janus colloid, measured at a distance r sp f, = 1.5R, 
from the center of the colloid. The diamonds correspond to an ideal- 
gas fluid with ap = 0, and the circles correspond to a DPD fluid with 
ap = 39 kgT/r c . The filled and the open symbols correspond to the 
cases with bounce-back and specular reflection BCs, respectively. 
Flere, the simulations are performed using a,A = am = 312 kpT/r c 
and the temperature T co ui is controlled directly at the colloid surface. 

distance from the center of the colloid, see Fig. [TO] The radial 
component of the velocity is measured at a particular angle, 
9, relative to the symmetry axis of the colloid. Within the 
statistical accuracy of our simulations, we observe a power- 
law dependence with r sp / 7 , which is consistent with an inverse 
cubic power, v rsph ~ 1 / r 3 ' h , for both an ideal-gas fluid with 
cif = 0, and a DPD fluid with non-zero ap. This behavior is 
unaffected by the different BCs used, and is in good agreement 
with the theoretical predictions-^ and MPC simulation s 24 : 39 . 
We also compute the radial velocity as a function of the angle 
0 at a distance r sp h = 1.5 R s from the center. The results are 
shown in Fig. HU for different fluid interaction strengths and 
different BCs. We find an asymmetric radial velocity distri¬ 
bution around the colloid from the cold cap pole (9 = 0°) to 
the hot cap pole (9 = 180°). This dependence is in qualitative 
agreement with the theoretical predictions of Ref. 21 . 

4 Discussion and conclusions 

The swimming velocity of a thermophoretic Janus colloid 
strongly depends on different fluid-colloid interactions and 
BCs. The choice of BCs (i.e., slip or no-slip) affects vis¬ 
cous friction exerted on the colloid by the fluid, while differ¬ 
ent strengths of fluid-colloid repulsive interactions alter both, 


the heat exchange between the fluid and the colloid surface 
and the near-wall density fluctuations illustrated in Fig. 0 
Thus, the repulsive strengths a ,4 and a,n also have an effect 
on the viscous friction for the case of bounce-back reflections. 
Perhaps, the simplest case is that with specular reflections in 
Fig. 0b), since viscous friction between the colloid and the 
fluid can be neglected. In this case, the colloid mobility in¬ 
creases with increasing interaction strength, or equivalently 
when fluid particles are pushed further away from the surface 
of the swimmer. We expect that a larger distance between fluid 
particles and the colloid surface should lead to a reduction of 
heat exchange between the colloid and the fluid. In Fig. 0b) 
for small a,/\ and am values, a fast heat exchange between the 
colloid and the fluid is expected. As the heat exchange is get¬ 
ting reduced for increasing fluid-colloid interaction strength, 
the swimming velocity is increasing. This indicates again that 
the temperature gradient between the cap and the fluid mainly 
determines colloid propulsion. 

The comparison of Figs. 0a) and 0b) for the stick and 
slip BCs have shown that the no-slip BCs also result in an 
enhancement of the swimming velocity for the low strengths 
of fluid-colloid interactions. This finding is rather counter¬ 
intuitive, since no-slip BCs lead to an additional friction on the 
colloid exerted by the fluid. Bounce-back reflections do not af¬ 
fect fluid particle distribution near the colloid in comparison to 
specular reflections, and therefore, it is plausible to expect no 
change in heat exchange between the colloid and the fluid at 
least through the heat conduction term in Eq. ©. The particle 
kinetic energy also remains conserved for both bounce-back 
and specular reflections. One should expect differences in the 
potential energy for the different collision rules. For exam¬ 
ple, bounce-back reflections may lead to a slight elevation of 
temperature (2 — 5%) near a wall in comparison to a specu¬ 
lar type of reflections, which has been found for the standard 
isothermal (non-energy-conserving) DPD^S,. A local increase 
of temperature near the colloid surface would reduce conduc¬ 
tivity between the colloid and the fluid, which would be con¬ 
sistent with an increase of the swimming velocity as discussed 
above for the fluid-colloid interactions. However, currently we 
cannot exclude that other effects are present and the interplay 
between viscous friction and heat exchange between the col¬ 
loid and the fluid for various parameters needs to be investi¬ 
gated in much more detail. 

Results for the ideal fluid have shown swimming trends 
qualitatively different from those for a non-ideal liquid. For 
instance, in case of bounce-back reflections the swimming ve¬ 
locity first decreases with increasing the fluid-colloid repulsive 
interaction strength for an ideal fluid (Fig. 0a)), while in the 
corresponding case for a non-ideal fluid (Fig. 0a)) the swim¬ 
ming velocity is increasing when the fluid-colloid interactions 
are getting stronger. In the case of specular reflections, the 
trends of an increase of the swimming velocity with increas- 










ing the fluid-colloid interaction strength (Figs. HJb) and[6lb)) 
are similar for both fluid types; however, for the ideal-fluid 
case the swimming velocity changes its sign, which means 
that the thermophoretic swimmer changes its swimming direc¬ 
tion. From Figs. [SJa) and[8jb) we can also conclude that the 
swimmer’s velocity in case of an ideal fluid becomes nearly 
independent of the repulsion strength and type of the fluid 
particle reflection, when a iA > 100 and am > 100, since such 
a repulsion strength is large enough to nearly push all fluid 
particles away from a layer of the cutoff radius r' c of the fluid- 
colloid repulsive interaction evidenced from fluid-density dis¬ 
tributions. Note that this occurs due to a high compressibility 
of an ideal fluid, while for a non-ideal fluid the layer of r' c 
still remains populated by fluid particles. Therefore, the heat 
exchange between the colloid and the fluid is expected to be 
affected by the fluid-colloid interactions much more for ideal 
fluids than that for non-ideal liquids. Another difference be¬ 
tween the ideal and non-ideal fluid cases is density gradients, 
which are much stronger for an ideal fluid than for a non-ideal 
liquid, even though temperature distributions are not drasti¬ 
cally different. For the case of a non-ideal fluid, the changes 
in density are within a few percent from an average fluid den¬ 
sity, while for the ideal-fluid case the density may change up 
to 30-40% from an average density. Thus, the fluid-density 
changes are more realistic in the non-ideal fluid case. How¬ 
ever, currently it is not clear how these differences between 
the ideal and non-ideal fluid cases lead to different swimming 
behavior of the thermophoretic colloid. 

A qualitative explanation for the behavior of a (homoge¬ 
neous) colloidal particle in a temperature gradient has been 
sketched in Ref.— for MPC simulations. A temperature gra¬ 
dient results in an inverse gradient of density and thus, the 
density around the cold side is larger than that at the hot side. 
Hence, a higher density on one side may result in a stronger 
interaction and lead to the colloid motion. A change in fluid- 
colloid interactions (e.g., repulsion or attraction) may invert 
this balance and force a colloid to move to an opposite di¬ 
rection. This idea is equivalent to having a pressure gradient 
across the colloid poles and the fluid-colloid interactions seem 
to provide a control for it. However, this argument for the gen¬ 
eration of a pressure difference is not completely conclusive, 
because for an ideal gas in local thermodynamic equilibrium, 
the pressure p = kgTp should be constant due to mechanic 
stability. 

Swimming of the Janus colloid toward the hot side in Fig. 
db) for a case of ideal fluid, specular reflections, and weak 
fluid-colloid interactions is consistent with this proposition. 
Note that in case of specular reflections no exchange of mo¬ 
mentum occurs between the fluid and colloid in the tangen¬ 
tial direction, and therefore, a driving force for the swimming 
colloid is likely to come from a pressure gradient across the 
colloid poles. Then, as we increase the repulsion between the 


fluid and colloid in Fig. |8jb), the pressure difference is turned 
around and the thermophoretic swimmer moves toward an op¬ 
posite direction. The comparison of the specular-reflection 
case in Fig. [8jb) to the bounce-back BCs in Fig. d a ) indi¬ 
cates that exchange of momentum between the fluid and col¬ 
loid in the tangential direction also contributes to the swim¬ 
mer propulsion, since the mobility of colloid is different for 
these two conditions. Note that no pressure differences are ex¬ 
pected between the specular and bounce-back cases, because 
both conditions lead to the same density distributions and the 
same exchange of momentum between the fluid and colloid in 
the normal direction. Currently, we cannot identify the swim¬ 
ming effect due to the tangential momentum exchange, but it 
is clearly present in these systems. 

The application of the idea above to the case of a non-ideal 
fluid is not so straightforward, because the density gradients 
in this case are much smaller than in the case of an ideal fluid, 
as already mentioned. In fact. Fig. [6JB) for specular BCs in¬ 
dicates that a colloid does not swim for the case of vanish¬ 
ing fluid-colloid interactions, which implies no pressure gra¬ 
dient across the colloid poles. As the repulsion between the 
fluid and colloid is increased in Fig. [6jb), the swimmer starts 
moving toward the cold side indicating that a pressure dif¬ 
ference across the colloid poles must have developed. The 
comparison of results in Fig. |6jb) and Fig. [6} a) for specular 
and bounce-back BCs, respectively, implies again that the ex¬ 
change of momentum between the fluid and colloid in the tan¬ 
gential direction must contribute to colloid’s swimming. Here, 
for the case of bounce-back BCs in Fig. [6|a), the colloid has 
a non-zero swimming velocity for vanishing fluid-colloid in¬ 
teractions. Thus, the effect of tangential momentum exchange 
on the propulsion of a thermophoretic swimmer needs to be 
investigated further. 

Finally, we would like to discuss the differences between 
our simulations with an ideal fluid and the MPC simulations 
of a similar Janus-colloid swimmer . 2 2 i — > 2 4 . Note that a direct 
comparison has not been intended. A seeming dissimilarity 
is the dependence of swimming velocity on the fluid-colloid 
interaction strength. In the present simulations the ther¬ 
mophoretic colloid swims toward the cold side if we increase 
the repulsive strength of fluid-colloid interactions, while in 
Refs. 20 ' 24 repulsive interactions between a colloid and a fluid 
result in the motion toward the hot side. A closer look at the 
details of the simulation setups reveals that the fluid-colloid in¬ 
teractions in these two studies may have a different meaning. 
In our setup, the repulsive interactions directly affect heat ex¬ 
change between the colloid and the fluid, since larger distances 
between them reduce the exchange of heat. In Refs. 20,24 the 
temperature in a thin layer near the colloid is controlled, and 
thus, the repulsive interactions affect the number of particles 
to be thermalized in this layer. Hence, in the present simu¬ 
lations repulsive interactions affect conductivity between the 





colloid and the fluid, while in Refs.— 24 such conductivity ef¬ 
fects are omitted, which is equivalent to a very high colloid- 
fluid conductivity such that a thin layer of fluid particles near 
the colloid receives heat instantly. Another difference between 
the two simulation setups is temperature control. In our simu¬ 
lations, generated heat is taken away at the cold side of the col¬ 
loid or far away from the colloid, while in Refs.^2 4 the excess 
heat is taken away uniformly from the whole fluid. This may 
result in different temperature distributions around the colloid 
affecting its swimming behavior. Finally, there exist a funda¬ 
mental difference between the simulation methods. In MPC, 
heat exchange and temperature gradients are sustained only 
through the kinetic energy of fluid particles, while in eDPD 
an internal energy is simulated explicitly. In fact, the internal 
energy is much larger than the contributions from kinetic and 
potential energies, since C v 1. Simulations with C v ~ 1 ap¬ 
pear not to be stable in eDPD, since then there is a chance that 
internal temperature of a particle may become negative, for 
instance, due to the random conductivity term in Eq. Q. The 
discussed reasons do not allow us to make a detailed compar¬ 
ison, which would require more consistent setups. 

In conclusion, we have presented simulations of the dynam¬ 
ics of a thermophoretic colloid for different fluid-colloid in¬ 
teractions and temperature controls. Different temperature- 
control strategies have a minor effect on the colloid swim¬ 
ming velocity. The fluid-colloid interactions have a strong ef¬ 
fect on the colloid behavior and directly affect heat exchange 
between the colloid surface and the fluid. Our results show 
that a reduction of the heat exchange leads to an enhancement 
of colloid’s thermophoretic mobility, since larger temperature 
gradients near the colloid surface are formed. The flow-field 
generated by the colloid appears to be dominated by a source- 
dipole contribution in agreement with the recent theoretical— 
and simulation 20 i 24 predictions. However, the differences in 
colloid’s mobility between the cases with non-ideal and ideal 
fluids and in comparison to the MPC simulations 20 24 are yet 
to be understood. We hope that this work will generate further 
efforts and discussions in this area of research. 
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